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Abstract In this study we numerically calculate the spatial 
profile of mechanical strain on self-assembled germanium 
(Ge) quantum dots (QDs) grown on a silicon (Si) substrate. 
Although the topic has been exhaustively studied, interest¬ 
ing features have not been explained or even mentioned in 
the literature yet. We studied the effect of the cap layer con¬ 
sidering two cases: capped QDs (where a Si cap is present 
above the Ge QDs) and uncapped QDs (where no Si is present 
above the Ge QDs). We observed that Ge in the capped QDs 
is more strained compared with the the uncapped QDs. This 
expected effect is attributed to the additional tension from 
the Si cap layer. However, the situation is opposite for the 
Si substrate, it is more strained in the uncapped QD because 
the Ge layer is less strained in this case. We also calculated 
the band-edge alignment for the electrons and holes. 

Keywords quantum dots • mechanical strain • uncapped 
quantum dot • capped quantum dot 

1 Introduction 

Mechanical strain plays an important role in electronic, opti¬ 
cal, and transport properties of semiconductor quantum dots 
(QDs) [1,2,3]. Epitaxial dots have a spatial strain profile 
caused by the mismatch of lattice parameters of the involved 
materials. This strain creates an elastic energy which is one 
of the determining factors in the formation of dots in the 
Stranski-Krastanov [4] mode. The strain also changes the 
potential profile experienced by charger carriers, necessary 
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to calculate the electronic eigenstates [5,6]. Several meth¬ 
ods to calculate the strain profile are available, which have 
been successfully applied to semiconductor nanostructures. 
The most used methods involve the continuous elasticity [7, 
8] or atomistic models [9, 10]. Several theoretical and exper¬ 
imental works on germanium (Ge)/silicon (Si) QDs can be 
found in the literature [1 1,12,13, 14], although none of them 
describes in detail the features of the strain and tension in 
uncapped QDs. 

In the present work, we studied the effect of the cap layer 
on the mechanical strain on self-assembled Ge QDs grown 
on silicon a Si substrate (Ge/Si QDs). We calculated the 
strain profile by numerically solving the equilibrium equa¬ 
tions in the continuous elasticity model using finite element 
method. The rest of this paper is organized as follows: in the 
next section, we present the model and describe the param¬ 
eters considered in the calculations. In Section 3, we show 
the results of our simulations. In Section 4, we present our 
conclusions. 


2 Theoretical model 

2.1 Linear elastic theory 

Axial symmetry is an approximation where a three-dimensional 
structure is generated by the revolution of a two-dimensional 
one. Dome (or lens-shaped) Ge QDs [15,16] belong to this 
case. This approximation requires that all physical param¬ 
eters of a crystalline lattice possess this symmetry, which 
means that all the relevant quantities of the problem do not 
depend on the angular coordinate (p of the cylindrical coor¬ 
dinates (r,(p,z). This approximation allows the problem to 
be solved in the two-dimensional plane r, z for r > 0 and it 
is very useful because it considerably decreases the compu¬ 
tational cost of the simulations. 
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The model used in this work is based on the numeri¬ 
cal solution of the equilibrium equations of the linear elastic 
theory [17,18]. This theory is suitable for treating materials 
with small strain and is based on the Hooke’s law, namely, 
F = —kx. Using a two-order tensor for strain Cij and stress 
Gij, Hooke’s law is expressed as [19]: 
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where Qj is the elastic constants matrix of a material with 
cubic symmetry. Conversely, if the material has a cylindri¬ 
cal symmetry, we need to perform a coordinate transforma¬ 
tion (equivalent to x r and y ^ (p) to obtain the Hooke’s 
law in cylindrical coordinates. In addition, by considering 
axial symmetry approximation, no quantity has angular de¬ 
pendence. Then, in cylindrical coordinates we have 
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where v is the displacement in the y—direction. Therefore, 
the Hooke’s law can be simplified as 
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The stress tensor is defined in terms of the external force 
as V • Gij-\-¥ = 0. Because we do not consider any displace¬ 
ment in the center of the mass, we obtain F = 0; therefore 


V-(T,-,=0. (4) 

Substituting Eq. (3) into Eq. (4) yields 

V^{Cij5eij)=0. (5) 

where the solution for the strain is represented as Scij. To 
obtain this solution, we numerically solve the above differ¬ 
ential equation [Eq. (5)] using the software Comsol Multi¬ 
physics. 



Fig. 1 Illustration of the dome-shaped geometry modeled as a spheri¬ 
cal cap. The dome is the top of a sphere whose radius is 43.8 nm. 


where £i is the components of the initial value of the strain 
needed to solve Eq. (5). The characteristic function x is de¬ 
fined as: 



1 if r G 
0 if r G ilsi 


( 6 ) 


where and Qsi refer to the space domains occupied by 
Ge and Si, respectively. The purpose of X is to take into 
account the assumption that the initial strain is zero in the 
Si domain and constant and nonzero in the Ge domain. In 
addition, the shear components of the initial strain are set to 
zero, i.e.. 


^rz — ^r(p — ^(pz — (7) 

Hereafter, we show how to calculate the initial values. In 
an epitaxial film containing QDs (with lattice parameter at) 
grown over a substrate (001) (with parameter am), the strain 
in the plane is biaxial and given by the mismatch in lattice 
parameter a, which is defined as [20] 

a = (8) 

Gi 

Where we used the the following assumption: 

£r = £cp = OJ. (9) 

In this condition, the stress tensor can be approximated by: 


(yr = (Jcp= 0. (10) 

Substituting Eqs. (9) and (10) into the Hooke’s law yields 

[ 20 ] 

e, = -2a^. (11) 

Gil 


2.2 Initial condition £ for the strain 

The components of the total strain are expressed as follows: 

e(p = X^(p^^e(p 
Cz = X^Z ^^z 

Crz = 


2.3 Lattice parameter and stress tensor 

The strain e^ir^z) gives the relative variation in the distance 
between nearest neighbors in the z—direction. In plane xy, 
the strain is the average of the other two components: 



(12) 
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Capped QD 





Fig. 2 Spatial profile of Cpir^z) [Eq. (12)] (on the left side) and 
(on the right side) for (a) Capped QD and (b) Uncapped QD. The strain 
is multiplied by 100 to be expressed in % and the color scale bar is the 
same for both graphics. 


The lattice parameter in plane xy is defined as [21]: 

aMz) = [\+er{r,z)] [Z«Ge + (l -Z)asi], (13) 

where uqq and are the lattice parameters of Ge and Si, 
respectively, and X is defined in Eq. (6). We also calculated 
the hydrostatic and biaxial components of the strain, which 
are defined as [20] : 

eh{r,z) =e^ + ep (14) 

and 

eb{r,z)=e^-ep, (15) 

respectively. The hydrostatic component reflects the varia¬ 

tion in the volume of the structure due to the strain. From 


Hooke’s law, the stress tensor is expressed as: 
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2.4 Band-edge alignment 

The difference in gap energy of the materials (in our case: 
Ge and Si) creates a band offset acting as potential well for 
the carriers. We represent this potential by Uc [for the elec¬ 
trons in the conduction band (CB)], and by Uy [for the holes 
in the valence band (VB)]. The hydrostatic component ek 
shifts this potential and the biaxial one brakes the degen¬ 
eracy of CB splitting the sixfold-degenerate A valleys into 
the fourfold A (4) and twofold A (2). In VB, the biaxial strain 
also splits the heavy and light hole bands. We write <Pa2 and 
0^4 for the energy levels of the electrons and and 
for the energy levels of the heavy and light holes [14]. Their 
dependency with the strain are given by 


^A 2 {r,Z,eh,eh) = Uc{r,z)+C^ 2 {eh,eb) (17) 

’^AA{r,z,eh,eb) = Uc{r,z)+CAA{eh-,eb) (18) 

'I'hh{r,Z,eh,eb) = Uy{r,z)+Vhh(eh,eb) (19) 

'l'ih{r,z,eh,eb) = Uy{r,z)+Vih{eh,eb) (20) 


The confinement potentials are defined as 
Uc{r,z) = Es+xUbc 
U v{r,z) = xUbv 

where = 1.17 eV is the energy gap of the Si at 0 K. The 
contributions of the strain in each potential are written as 
[22,23,24]: 

2 

~ “ 1 “ 

CAA{eh,eb) = ahCh - -de^ 

^hhi^hj^b) ~ ^^b 

yih{eh,eb) ^ + 2 ^^^ + 2 ~ 9beb)^ + ^Aobet, 

where ah = ac — a^. The importance of these band-edge en¬ 
ergies is that they enter as input in the Schroedinger equa¬ 
tion in order to calculate the wave-functions of the carriers. 
Therefore, a minimum in the potential can represent a local¬ 
ization of the wave function. 


2.5 Calculation parameters 

With regard to the geometry of the QD, we relied on the 
results obtained by Magalhes-Paniago et. al [25]. They ob¬ 
served a dome-shaped geometry with average radius R = 
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Quantity 

Si 

Ge 

a (A) 

5.431 

5.657 

Cii (10“ Pa) 

1.66 

1.29 

Ci 2 (10“ Pa) 

0.64 

0.48 

C 44 (10“ Pa) 

0.8 

0.67 

(eV) 

-3.133 

-1,786 

(eV) 

-3.1 

-4.54 

d(eV) 

8.6 

9.4 

4o (eV) 

0.044 

0.296 

b(eV) 

-2.1 

-2,86 


Table 1 Constants and parameters of Si and Ge. The values were ob¬ 
tained from Ref. [12]. 


32.3 nm and height H = 14.2 nm (Fig. 1). We also consid¬ 
ered a quantum well (also called wetting layer) with thick¬ 
ness L = 1 nm between z = 0 and z = L whereas the QD is 
located between L < z < L-\-H. Using nanometer (nm) as 
the units for space, the Si substrate comprises the domains 
—500 < z < 500 and 0 < r < 500. The wetting layer com¬ 
prises 0 < r < 500 and 0 < z < 1. The calculation was done 
only in the region r > 0. From Eq. (8) and Table 1 we have 

a = = _o.039. 

ClGe 


3 Results 

Figure 2(a) shows the spatial profile of the strain in plane 
ep{r^z) on the color chart for the capped QD. A negative 
strain indicates that the distance between the atoms has de¬ 
creased. Because Ge has a lattice parameter that is greater 
than that of Si, Table 1 , the epitaxially grown Ge layers have 
a negative strain in plane vy of the order of approximately 
—4%. Just as the Si layer compresses the Ge QD, the Ge 
QD expands the Si layer, which is what we observed in the 
Si layers immediately above and below the Ge QDs where 
the strain reaches positive values higher than 0.5%. 

Figure 2(b) shows the spatial profile of ep{r^z) for the 
uncapped QD. The analysis follows the same trend: nega¬ 
tive strain at the QD and positive one for the Si layer be¬ 
low the dot. However, the difference is that no Si layer is 
presented above the Ge QD, which results in a remarkable 
distinction: the strain in Ge is lower in module, i.e., while 
at the capped QD the strain is approximately —3% at half 
height, it is approximately —2% at the uncapped QD. We 
can understand this by considering that in the capped case, 
both the substrate and Si cover compress the Ge. In the un¬ 
capped case, only the substrate compresses the QD, and at 
the superior interface, the Ge is free; thus, its original lattice 
parameter tends to be reached more quickly. The behavior 
is different in the Si substrate, which is more strained un¬ 
der the uncapped dot. This condition can be understood as a 
reaction of the Ge dot, which in turn is less strained and pos¬ 
sesses a more different lattice parameter that creates a larger 




r (nm) 


Fig. 3 Spatial profile (in color) of the resultant lattice parameter in the 
plane ar{r,z) [Eq. (13)]. The blue lines indicate constant values of ar. 
(a) Capped QD. (b) Uncapped QD. The color scale bar is the same for 
both graphics. 


Strain on the Si layer. In the capped case, the Ge dot is more 
strained, which makes its lattice parameter more similar to 
the Si layer, creating a small strain on the Si layer. 

The z— component of the initial strain is proportional to 
the component in the plane but with an opposite sign. When 
the first layer of Ge is bound to the Si layer, compression 
exists in the plane of Ge, but no tension is present in the 
z— direction perpendicular to the plane. Therefore, the com¬ 
pression in the plane (ep < 0) implies an expansion in the 
z— direction (e^ > 0) as if Ge tries to keep the volume of 
the unit cell costant. Figure 2(a) shows the behavior of 
for the capped QD. One can see that a small compression 
(e^ ^ —0.5 %) exists in the Si layer under the QD and a large 
expansion (e^ ^2%) exists inside the QD in response to the 
compression in the plane. In the same manner, the Si cap 
layer also exhibits compression (^ — 1 %). In the uncapped 
case, shown in Fig. 2(b), no Si cap layer exists to exert an 
additional stress on the Ge QD, which makes the uncapped 
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Fig. 4 Spatial profile (in color) of the tension in the plane Op(r,z) [Eq. (21)]. (a) Capped QD. (b) Uncapped QD. The color scale bar is the same 
for both graphics. 
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Fig. 5 Band-edge alignment in the z—axis (r = 0). Solid lines refer to 
the capped QD and dot-dashed lines represent the uncapped QD. The 
silicon substrate is the region z < 0, the wetting layer is 0 < z < 1, the 
QD is 1 < z < 15. (a) Conduction band (CB). Blue lines are <Pa 2 [Eq. 
(17)] and red lines are <Pa 4 [Eq. (18)]. Black dashed line is Uc. The cap 
layer is z > 15. (b) Valence band (VB). Red lines are Whh [Eq. (19)] and 
blue lines are Wih [Eq. (20)]. Black dashed line is t/y. 


QD less strained than the capped QD, especially near its top 
surface. This condition makes the uncapped QD have a sim¬ 
ilar lattice parameter to its bulk value. Consequently, the un¬ 
capped QD exerts more tension on the Si substrate compared 
with the capped case. 

The lattice parameter in the plane [Eq. (13)] is shown 
in Fig. 3(a). As expected, the lattice parameter in the Ge is 
smaller than its bulk value. At the top of the dot, ay reaches 
the maximum value for Ge, which is 5.5 A. In the same 
manner, in the Si above and below the dot, the parameter 
is larger than its bulk value (5.431 A). In both regions, Uy 
reaches values greater than 5.45 A. Figure 3(b) shows the 
profile of ay for the uncapped dot. The lattice parameter in 
Ge attains higher values in the uncapped dot than that in the 
capped one. For instance, at the top of the uncapped dot, 
Uy attains values exceeding 5.6 A. However, the maximum 
value of the capped dot (5.5 A) is only the minimum value 
of the uncapped dot. The Si also shows a difference. In the 
uncapped dot, the Si below the dot attains higher parameters 
of the order of 5.5 A. 

Figure 4(a) presents the calculated stress tension in the 
plane, 

defined as analogous to ep. The components Gy, G(p, and G^ 
are defined in Eq. (16). Figure 4(b) shows Gp{r^z) for the 
uncapped dot. In the capped dot, the compressive tension in 
Ge is greater than that in the uncapped one. Conversely, in 
Si, the tension of traction (which tends to expand) is greater 
at the uncapped dot than that in the capped one. This is ex¬ 
pected whereas in the capped dot, Ge undergoes more com¬ 
pression and thus manages to exert less traction (or tension 
of expansion) in Si at the capped QD. 
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Figures 5(a) and 5(b) show Uc and Uy (black dashed line) 
without strain contribution. In this situation, the electron is 
localized in the Si layers (substrate or cap layer) while the 
holes are confined in the QD. Figure 5(a) shows the graphs 
of the electron energy levels A (2) and A (4) in the CB along 
the z—axis. For the capped QD, the minimum in the Si layer 
is in the A (2) band while inside the QD, the minimum is in 
the A (4) band. The split of these two bands are ~ 83 meV 
in the substrate and ^480 meV inside the QD. In the un¬ 
capped case, the minima in the substrate and inside the QD 
remain unchanged although the energy split inside the QD 
decreases to ^ 336 meV. Additionally, the band edge inside 
the uncapped QD has its minimum in the base of it, near the 
substrate. Other interesting feature is that the A (4) minimum 
inside the uncapped QD is ^a 4 (z = 2) — fpAiiz = — 1) = 83 
meV higher than the A (2) minimum in the substrate. In the 
capped case, this difference is only 14 meV. This means 
that the electron wave-function (which is localized in the Si 
layer) can penetrate much less inside the uncapped QD than 
in the capped one. A smaller penetration implies a lower su¬ 
perposition with the hole wave-function. This can be impor¬ 
tant to understand the different behavior of the band align¬ 
ment (type I or type II) already studied for this system [14, 
26,27]. 

Figure 5(b) shows the energy levels of the heavy and 
light hole in the VB. An interesting feature is the energy split 
of the bands, Wih — ^hh at z = 7.5 nm, which is 384 and 240 
meV for the capped and uncapped cases, respectively. As 
the two holes are already confined within the QD, the effect 
of the strain does not change considerably their eigenstates. 
However, in experiments where an external strain is applied, 
the two bands (hh and Ih) shifts with different speeds. In 
the uncapped case, an energy split smaller than the above 
values can cause an anti-crossing of the two bands [28] with 
a smaller external strain. 

4 Conclusions 

We have calculated the spatial profile of mechanical strain 
in self-assembled Ge QDs grown on Si matrix. The Ge QD 
experiences a compressive biaxial strain in the vy—plane, 
which in turn creates a tensile strain on the z—plane. This 
effect is less intense in the uncapped QD because only one 
Si layer (substrate) compresses the Ge QD. Conversely, the 
Si substrate and cap layer experiences tensile strain on the 
vy—plane and compressive strain on the z—plane. In addi¬ 
tion, the Si substrate is more strained in the uncapped case 
because the Ge uncapped dot is less strained. All these ef¬ 
fects are reflected on the final lattice parameter and the stress 
tensor. We also calculated the band-edge alignment of the 
conduction and valence bands. The strain induced poten¬ 
tials suggest a smaller wave-function superposition in the 
uncapped QD when compared with the capped one. 


References 

1. J. Stangl, V. Hoi, G. Bauer, Review Modern Physics 76 (2004) 725- 
783. 

2. F. Boxberg and J. Tulkki, Reports on Progess in Physics 70 (2007) 
1425-1471. 

3. R. S. Jacobsen, K. N. Andersen, P 1. Borel, J. Fage-Pedersen, L. H. 
Frandsen, O. Hansen, Nature 441 (2006) 199-202. 

4. G. Medeiros-Ribeiro e R. S. Williams, NanoLetters 7 (2007) 223. 

5. M. Grundmann, O. Stier, and D. Bimberg, Physical Review B 52 
(1995) 11969-11981. 

6. A. D. Andreev and E. P O’Reilly, Physical Review B 62 (2000) 
15851-15870. 

7. B. Jogai, Journal of Applied Physics 90 (2001) 699-704. 

8. A. Cros, N. Garro, J. M. Llorens, A. Garca-Cristbal, A. Cantarero, 
N. Gogneau, E. Monroy and B. Daudin, Physical Review B 73 (2006) 
115313. 

9. O. Stier, M. Grundmann, D. Bimberg, Physical Review B 59 (1999) 
5688-5700. 

10. L. He, G. Bester e A. Zunger, Physical Review B 70 (2004) 
235316. 

11. A. 1. Yakimov, A. A. Bloshkin and A. V. Dvurechenskii, Physical 
Review B 81 (2010) 115434. 

12. K. Brunner, Reports on Progress in Physics 65 (2002) 27. 

13. M. El Kurdi, S. Sauvage, G. Eishman, and P Boucaud, Physical 
Review B 73 (2006) 195327. 

14. P E. Gomes, E. likawa, E. Cerdeira, M. Larsson, A. Elfving, G. 
V. Hansson, W.-X. Ni, and P-O. Holtz, Applied Physics Letters 91 
(2007) 051917. 

15. E. M. Ross, R. M. Tromp, M. C. Reuter, Science 286 (1999) 1931. 

16. G. Medeiros-Ribeiro, A. M. Bratkovski, T. I. Kamins, D. A. A. 
Ohlberg, R. S. Williams, Science 279 (1998) 353. 

17. L. D. Landau e E. M. Lifshitz, Theory of Elasticity, 3rd Edition, 
Landau and Lifshitz Series Course of Theoretical Physics, volume 7, 
Elsevier, (1986). 

18. Martin H. Sad, Elasticity: Theory, Applications an Numerics, El¬ 
sevier Academic Press (2004). 

19. L. A. Montoro, G. Medeiros-Ribeiro e A. J. Ramirez, Journal of 
Physical Chemistry C 114 (2010) 12409-12415. 

20. M. K. Kuon T. R. Lin and K. B. Hong, Journal of Applied Physics 
103 (2008)073705. 

21. D. Jalabert, J. Coraux, H. Renevier, B. Daudin, M.-H. Cho, K. B. 
Chung, D. W. Moon, J. M. Llorens, N. Garro, A. Cros, and A. Garca- 
Cristbal, Physical Review B 72 (2005) 115301. 

22. C. Pryor, J. Kim, L. W. Wang, A. J. Williamson, and A. Zunger, 
Journal of Applied Physics 83 (1998) 2548. 

23. O. G. Schmidt and K. Eberl, Y. Rau, Physical Review B 62 (2000) 
16715. 

24. M. Califano and P. Harrison, Journal of Applied Physics 91 (2002) 
389. 

25. R. Magalhes-Paniago, G. Medeiros-Ribeiro, A. Malachias, S. Ky- 
cia, T. I. Kamins, and R. Stan Williams, Physical Review B 66 (2002) 
245312. 

26. M. Larsson, A. Elfving, W.-X. Ni, G. V. Hansson, P. O. Holtz, 
Physical Review B 73 (2006) 195319. 

27. S. V. Kondratenko, A. S. Nikolenko, O. V. Vakulenko, M. Ya 
Valakh, V. O. Yukhymchuk, A. V. Dvurechenskii and A. I. Nikiforov, 
Nanotechnology 19 (2008) 145703. 

28. P. E. Gomes, M. P. E. Godoy, M. K. K. Nakaema, E. likawa, T. E. 
Lamas, A. A. Quivy, and J. A. Brum, Physica Status Solidi (c) 1 No. 
3 (2004) 547. 



